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Abstract 

This paper introduces a new model of continuous opinion dynamics with 
random noise. The model belongs to the broad class of so called bounded 
confidence models. It differs from other popular bounded confidence models 
by the update rule, since it is intended to describe how the single person can 
influence at the same time a group of several listeners. Moreover, opinion 
noise is introduced to the model. Due to this noise, in some specific cases, 
spontaneous transitions between two states with a different number of large 
opinion clusters occur. Detailed analysis of these transitions is provided, with 
MC simulations and ME numerical integration analysis. 

1. Introduction 

Models of opinion dynamics are very popular in modern sociophysics (see 
recent reviews 

BBS)- 

An interesting subset of them are the models with 
Bounded Confidence (BC models) 0, 0, 0, 3- 

In these models, the opinion 
exchange takes place only when the difference between two opinions is below 
the confidence bound (also called as threshold or tolerance). This is a rea- 
sonable consideration because if the minds of two people are very different, 
it is difficult for them to convince each other of something. Sometimes it is 
even hard for them to talk to each other. As a result of opinion exchange, 
one agent can change the opinion of another agent, or they can convince 
each other. There are many different types of opinion exchanges, including 
exchanges between more than two agents. BC models are commonly used 
to simulate the evolution of opinion distribution in a set of agents. Depend- 
ing on tolerance, the simulation results can give consensus (one big cluster 
in an opinion space) for big Tolerance, polarization (two big clusters) for 
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smaller Tolerance, or clusterization (three or more clusters) for even smaller 
Tolerance. 

Although several extensions based on continuous opinions and bounded 
confidence have been proposed and analyzed little attention has been paid 
to the instability of opinion caused by influences from outside 12j, |l3|, |l4 



Hence, the simulation in most BC models ends when static configuration is 
reached (consensus, polarization, or clusterization) and nothing else happens. 

It is noticeable that even in countries with stable democracy we can ob- 
serve oscillations in opinion distribution, some of them are very strong. This 
result cannot be obtained in simulations if only the interactions between 
agents is taken into account. The model described in this paper tries to 
incorporate influence from outside. 

I understand that many external factors can influence an agent's opinion, 
such as dramatic event, Mass- media or others. This influence generates some 
unpredictable opinion changes. It may also be regarded as free will opposing 
the conformist character of opinion exchange. I assumed that the opinion 
of such an influenced agent may change to a completely different one, as in 



models with discrete opinions where it is known as " contrarians" [8|, [9|, [10|, 111 



To simulate these various unpredictable changes noise was added to that 



model [12|, [13|, As in [13l . Il4| I decided to change from time to time, 



an opinion of one randomly chosen agent to a new opinion, chosen from 
uniform random distribution between and 1. Probability of this change 

Another type of noise 
but this did not 



16, 17 



will be described by noise parameter p [15l 
was introduced in BC models with continuous opinions 18 
significantly affect the simulation results. 

In the next section, a new BC model with continuous opinions and random 
changes in an agent's opinion (noise parameter p) will be introduced. 



2. Model 

Consider a set of N agents. Each agent is connected with the others (such 
a structure can be described by a complete graph) and has its own opinion 
which is represented by a real number between and 1. 

1. Randomly choose one agent from set A = {a\, a n } denoting it's 
opinion by S*. 

2. Randomly choose L agents from the rest of the set. Their opinions will 
be subset {Si} of A where i — {1, L}. These agents will be listeners. 
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3. For each ie{l,...,L} if \St - S*\ <T ^ S' t = \{S* + Si) . 

4. With probability p, randomly choose one and only one agent from set 
A and change its opinion value to a new randomly chosen value from 
division (0, 1). 

5. Back to 1. 

As described, the parameters of the simulation are: L - number of listen- 
ers, p - noise parameter, T - tolerance and N - number of agents. 

The model proposed here is quite similar to the most popular BC models, 
i.e., the Deffuant (D) model [Ej and that of Hegselmann and Krause (HK) 
Opinions take real values in an interval [0, 1] and each agent, with opinion 
S*, interacts with agents whose opinions lie in the range [S* — T, S*+T\. The 
difference is given by the update rule: chosen agent does not interact with 
one of its neighbors, like in D model, but with L compatible neighbors at 
once, similarly to HK model. However, on contrary to HK model, opinions 
of L neighbors are changed, instead of changing an opinion of chosen agent. 
This means that one agent influences simultaneously L compatible neighbors, 
instead of being influenced by them. Differences between these three models 
can be viewed shortly in the following way (see also [H): 

• Deffuants model describes the opinion dynamics of large populations, 
where people meet in small groups, like pairs. 

• HK model describes formal meetings, where there is an effective inter- 
action involving many people at the same time. 

• The model proposed in this paper is intended to describe how the single 
person can influence at the same time (during a formal meeting) a group 
of several listeners. 

As mentioned in the Introduction, to describe the various unpredictable 
changes, the noise was introduced to the model just like Pineda et al did 
with the Deffuant model. 

It's very important that results obtained with Deffuant model with noise 
are the same as with this model for L = 2. However there's slight difference 
between update rules but it has no impact on the MC results, and has nothing 
to do with analytical approach. Hence model described in this paper could 
be treated as generalisation of Deffuant model with number of interacting at 
one time persons as a parameter. 
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Due to the noise, the system never reaches the final fixed point, but rather 
dynamic equilibrium. Moreover, after some time the opinion distribution is 
independent of the initial conditions unlike in noiseless BC models case - if 
there are any two different initial distributions, for example (Q (e.q. uniform) 
and R (e.q. normal)), and p Q = p R > 0, L Q = L R , T Q = T R , N Q = N R , after 
some number of steps we cannot distinguish between these two systems. 
Surely their distributions in any given moment will be different, but if a 
certain timespan is given, their statistical properties will be the same. All 
the simulations were made after the system reached dynamic equilibrium. 

3. Results 

3.1. Monte Carlo simulations 

n(O) andae hp0| - Y 
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Figure 1: Opinion distribution n{0) could be approximated by / (O) w ae\~P°\~' 1 where 
a, (3 and 7 are the factors. 



• IfL>0,p = 0, and T — 1, the consensus is reached very quickly. 

• If L = 0, p > 0, and Te < 0, 1 >, there is no information exchange and 
uniform random distribution of opinions appears. 

• If L > 0, p > 0, and T — 1, the set of agents does not reach consensus 
but rather stays in dynamic equilibrium with one big opinion cluster. 
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Figure 2: Different shapes of opinion distribution n(0) for different For smaller 
maximum is higher and standard deviation of distribution is greater than for greater (p 
- noise parameter and L - listeners number). 



SD(p/L) 




Figure 3: Standard deviation of opinion distribution SD(j-) for N = 1600 and T = 1. 



For T — 1, distribution of the opinions can be approximated by: 

/ (0) « ae'-^l-r (1) 

rule, where a, (3 and 7 are the factors that depend on model parameters 
(see fig. [1]). For T = 1 and iV = const., standard deviation (SD) of opinions 
depends only on £ - the greater is factor £ the greater is SD (see figs. [2] and 
fig. [3]). As it can be seen, regardless of the force attracting the agents to the 
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Figure 4: Number of clusters in opinion space (in other words, number of modes in opinion 
distribution) for several values of the tolerance factor T, L = 16, p = 0.16, N = 1600. 



center, there are still some agents spread on the whole opinion space. 

For T < 1, as in other BC models, opinion fragmentation occurs. The 
smaller T is, the greater number of clusters occur (see fig. H|). Assuming 
X = T^f, the number of large opinion clusters C corresponds to 

C ps [AX + B] (2) 

rule (see fig. ED, which is more accurate than 



C « [X] (3) 

proposed by Deffuant [5]. Using X rather than T where X = is 
more suitable for the presentation of simulation results, so it is more often 
used in this paper. Most importantly the number of clusters C changes 
continuously, not discretely but there are evident steps. There are regions 
where C is far from the integer. These regions in the space of X are most 
interesting because criticality appears there, according to bifurcation points 
in 14(. Due to instability in these regions, spontaneous transitions between 



states with different numbers of big opinion clusters appear. 



6 



o 



C(X) and [AX + B] 



J 



1 



C(X) 
[AX + B] 



Figure 5: Number of opinion clusters for L = 64, p = 0.64, N = 1600, 10 3 MCS, 10 
simulations, as a function of inverse opinion threshold X = ^ and C « [AX + £] ap- 
proximation. 
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Figure 6: Number of large opinion clusters C(X) for different number of agents, only first 
unstable region, L = 16, p = 0.16, 10 4 A/CS\ 10 simulations. 



The Deffuant model for continuous opinion dynamics under the presence 
of noise has been studied recently [14]. They were able to derive a master 
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Figure 7: |1 >^ |2 > transition L = 16, p = 0.16, X = 1.66, N = 1600. Points a, b, c, d, 
e, f on upper panel of figure correspond to the six bottom panels, respectively. 



equation for the probability density function which determines the individ- 
uals density or distribution in the opinion space. Moreover, they have also 
found that in the noisy case the asymptotic steady-state probability distribu- 
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Figure 8: |2 >^ |1 > transition L = 16, p = 0.16, X = 1.66, N = 1600. Points a, b, c, d, 
e, f on upper panel of figure correspond to the six bottom panels, respectively. 



tions reached by Monte-Carlo simulations might not coincide with the ones 
obtained from the master equation 14(. This takes place for finite systems 
and is caused by perturbation introduced by noise. 
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Figure 9: Fragment of typical simulation outcome |2 >— > |1 > transition in points a, b 
and |1 >-> |2 > transition in points c, d L = 16, p = 0.16, X = 1.66, N = 1600. 
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Figure 10: Number of transitions t(X) per 10 3 MCS as a function of inverse opinion 
threshold X = ^ L = 64, p = 0.64, AT = 1600, 10 3 MCS 10 simulations. 



Observed deviations were more pronounced in the case of being close to 
a bifurcation point. 

In this paper we study a new model, which differs slightly from the Def- 
fuant model, yet belongs to the class of bounded continence models. It occurs 
that different BC models without noise exhibit very similar behavior [l|, 0]. 
Therefore, one could expect similar behavior under the presence of noise. In- 
deed results obtained using Monte Carlo simulations for D model with noise 



14| agrees with results obtained in this paper. Moreover, due to the simi- 
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Figure 11: Number of transitions per 10 3 MCS, t(X) for different number of listeners L, 
only first unstable region, p = 0.08, N = 1600, 10 4 MCS 10 simulations. 




Figure 12: Number of transitions per 10 3 MCS, t(X) for different noise parameter p, only 
first unstable region, L = 8, N = 1600, 10 4 MCS 10 simulations. 



larities between BC models, results obtained here suggest that spontaneous 
transitions appearing in a bifurcation point might be responsible for the in- 
consistency between analytical results and simulations observed in [H, 0]. 
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t(X) transitions per 10 3 MCS for different N 
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Figure 13: Number of transitions per 10 3 MCS, t(X) for different number of agents, only 
first unstable region, L — 16, p — 0.16, \Q A MGS 10 simulations. 

Let us now describe steady states of our system more carefully. For the 
description of the system's state, \k > notation will be used, where k denotes 
the number of big opinion clusters. 

• For Xe < 1, Xi — A >, the system has one big opinion cluster which 
is in |1 > state, where A = A(N,L,p) is a monotonically decreasing 
function of the total number of agents N. For N — > oo A — > 0, which 
is usual behavior in the critical point (see fig. E]). 

• For Xe < X\ — A,Xi + A > (first unstable region, see fig. E]), spon- 
taneous transitions between one and two big clusters occur (See fig. [7J 
fig. M), fig. EE It can be denoted as |1 >— > \2 > for one to two 
clusters transitions and |2 >— > |1 > for opposite. Finally there is 
|1 >— 7- |2 >— > |1 > cycle, where time intervals between transitions are 
unpredictable. 

• For Xe < X± + A, X 2 — A > the system has two big clusters, where 
the second critical point X = X 2 . 

• For Xe < X 2 — A,X 2 + A > (second unstable region), spontaneous 
transitions between two and three big clusters are observed. |2 >— > 
1 3 >— > 1 2 > cycle occur. 
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Figure 14: Distribution of transitions for parameters L = 4, p — 0.04, N = 1600. Results 
for only one simulation are presented and in this case observation time is 10 7 MCS. Disc 
line is T| 1> distribution where T|i> is 1 > state lifetime or time between |2 >— > |1 > and 
1 >^ |2 > transitions, circle line is T| 2 > distribution where T| 2 > is 2 > state lifetime or 
time between |1 >— > |2 > and |2 >^ |1 > transitions. 



• Generally, for Xe < X^-i + A, Xj, — A > opinions are fragmented into 
k clusters (k- modal opinion distribution). 

• For Xe < X^ — A,X& + A >, spontaneous transitions between k and 
k + 1 big clusters are observed and \k >— > \k + l >— > \k > cycle appears. 

It is surprising that such a simple model can simulate such a complex 
behavior. Once again it should be mentioned that spontaneous transitions 
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occur only in critical regions around the bifurcation points Xe < X^—A, + 
A > and that for N — > oo A — > 0, which is usual in the critical point. 
However, it should be noticed that in social systems N can rarely be treated 
as infinite and thus A > 0. Within the proposed model, 'real life' takes place 
in the critical region. 

Let me now examine the spontaneous transition's mechanism. 

• For |1 > state in the first unstable region there is one big cluster in 
the center and two small clusters near and 1 (see panel denoted by 
(a) in fig. [7]). The position of the central cluster is about 0.5 but it 
oscillates very strongly, and the oscillations are larger for smaller N 
and are vanishing for iV — > oo. Sometimes, when it goes far to one 
of the sides, the opposite small cluster grows very fast and becomes a 
second big cluster immediately and |1 >— > |2 > transition takes place 
(see fig. [7J and fig. [9]). 

• |2 >— > |1 > transition (See fig. [SJ, fig. E]) is different and more rapid. If 
the set is in |2 > state and it is in the first unstable region, there are 
two large clusters whose positions oscillate slightly. Sometimes they get 
so close to each other that their tails begin to interact and attract each 
other. As a result they become closer and closer, and finally joining 
into one big cluster. 

• In the next unstable regions the mechanism is similar. In \k >— > \k + 
1 > transition, a new cluster is created between two other clusters. In 
\k+ >— > \k > transition, two adjacent clusters join into one. 

Although the exact moment of spontaneous transition cannot be pre- 
dicted, its frequency t (average number of transitions per 10 3 MCS) can be 
measured. There are several maximums of transition frequency, exactly in 
the centres of unstable regions (float part of C(X) fh 0.5) (See fig. [TU]). This 
is logical because the instability of an opinion's distribution is greatest in 
such places, so even a small perturbation can cause the spontaneous reor- 
ganization of opinions in the whole set. In each unstable region the shape 
of t(X) can be alternately approximated by Gaussian distribution, and also 
follows the rule: 
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where \i = p(L, p, k, N). It can be seen (fig. [TD]) that /x is a decreasing 
function of k. This is understandable due to the fact that for a greater X 
(greater k and smaller T) an agent can interact with less number of listeners 
than in the case of a smaller X (smaller k and greater T), because they 
are out of confidence bound. One can say that effective L is smaller and 
analogous to fig. [TTJ the smaller L the fewer transitions. 

As mentioned above, the average number of transitions per time unit de- 
pends not only on X (or T), but also on parameters p, L and N, t(L, p, X, N). 
It has been already shown that factor £ determines the shape of n(0). On 
the other hand, for larger p or L, the transitions occur much more often (see 
fig. [T2] ITT]) . It is easy to understand why. For greater p, & is also greater, 
hence the clusters are wider and it is easier for them to interact. When L 
is greater, the fluctuations are also bigger. This is because when one agent 
whose opinion is quite rare speaks to many other agents, it can convince 
many agents to its rare opinion. As a result this rare opinion gets stronger 
and begins to attract many other agents. So the fluctuation grows and can 
cause the transition very easily. It's also worth to notice that for greater p 
critical region shifts towards greater X and for greater L it shifts in opposite 
direction. This phenomena occurs due to changes in shape of n(0) and for 
constant 4 there is no shift. Of course it should be mentioned here once 
again that for L = 2 this model behaves identically to Deffuant model with 
noise. 

Of course, as mentioned above, spontaneous transitions occur only in crit- 
ical regions Xe < Xk — A, X^ + A > and A = A(N, L, p) is a monotonically 
decreasing function of the system size N (see fig. [6]). Such behavior is typical 
for critical phase transitions. Dependence on N is also visible in t(X) (see 
fig. [TBI . The more agents present, the less transitions occur. Again, this can 
be easily understood. With more agents, the distribution of opinion is more 
stable because it is hard to obtain such big fluctuations as with a small num- 
ber of agents. Although for an infinite set there will not be any transitions 
identically as in [14| ,but social systems are finite and transitions may occur. 

An analysis of the distribution of time between transitions (i.e state's 
lifetime) r (see fig. [T4"j) gives more interesting details. The parameters of 
the system were set to L = 4, p = 0.04, N = 1600 and the observation 
time was 10 7 MCS with 1 simulation. Simulations were provided for several 
different X values in the first unstable region. The results are presented for 
X = 1.61, 1.63, 1.64, 1.65. It is clear that each state has its own characteristic 
lifetime scale - there are maximums on the lifetime histograms. The position 
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of the maximum depends on X. Where probabilities of two states are similar, 
the maximums of t([log 2 (T)]) distribution for higher and lower state are also 
close to each other (See fig. [TU b, c). When these probabilities differ, the 
distributions are also further and characteristic lifetimes differs more See 
fig. [Ha, d. 

3.2. Master Equation results 

Master equation for this model when L — 1 is: 



= (l-p)[ 2 [ dO'P{20-0\t)P{0',t) 

ot J\0-0'\<T/2 

P(0,t) f dO'P(0',t)]+p[P a (0)-P(0,t)}, (5) 

J\0-0'\<T 

and is almost identical to that one for Deffuant system with noise derived 



by Pineda |14j. In fact to get that equation we need to divide first part of 
Pineda's equation by two. This is due to fact that in Deffuant model two 
agents are changing their opinions at one step but in this model for L = 1 
there's only one. 

To figure out behavior of this equation I numerically integrated it using 
fourth order Runge-Kutta method. From analysis of stability comes out 
that for some values of X solution of this equation is unstable. I always 
started from uniform opinion distribution P(0, 0) = 1 and after reaching 
an asymptotic solution I introduced perturbation gP p (0), where P'(0,t) = 
(1 — g)P(0,t) + gP p {0) was the perturbed distribution. There were two 
types of perturbation: symmetric P s {0) = 2(1 — \20 — 1|) and asymmetric 
P a (0) = 20, both of them are normalized. It occurred that |1 > state 
is immune to symmetric perturbation and needs asymmetric one to make 
|1 >— > \2 > transition. In |2 > state case situation is opposite and symmetric 
perturbation is needed to make |2 >— > |1 > transition. 

Therefore, in the case of |1 > I was introducing asymmetric perturbation 
P a and I was checking how big g has to be, to make |1 >— > \2 > transition. 
If the solution was |2 > I was introducing symmetric perturbation P s and 
I was checking how big g has to be, to make |2 >— > |1 > transition. The 
dependence between the amount of perturbation p needed for transition and 
the parameter X is presented in FigJT5l 

It can be seen that for some values of X stability of solution is very low 
and it can switch to another solution very easily but opposite switch requires 



16 



perturbation p needed to make transition 
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Figure 15: Amount of perturbation p needed to make transition 



much greater perturbation. However, there are values of X where this switch 
requires perturbation of the same strength (the same p) in both directions, 
but still have to remember that in direction |1 >— > |2 > we need asymmetric 
perturbation, and for |2 >— > |1 > symmetric one. Because those two types 
of transitions are equally easy in points of equal p their should be analogous 
to the maximums of t(X) from MC simulations, and as can be seen in figJT5l 
they are very close. 

4. Discussion 

The model proposed in this paper is a simple bounded confidence (BC) 
model on a complete graph with noise p added to simulate outside influences 
as well as the free will and unpredictability of individual agents. One ran- 
domly chosen agent can communicate with L other randomly chosen agents 
according to the BC rule. Then, with probability p, the opinion of one ran- 
domly chosen agent is changed to a random number between and 1. As 
in the other BC models, clusterization of opinions occurs. The new qual- 
ity that appears due to the introduced noise is the presence of spontaneous 
transitions between different numbers of clusters. They take place for some 
specific values of X = ^ where T is tolerance. 

As is usual in the case of critical phenomena, in the proposed model 



|2|->|1> symmetric 
|1 1->|2> asymmetric 
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competition between two opposite forces is present. Due to competition 
between these forces, the system is in dynamic equilibrium. However, in 
some cases (for specific T values, close to the bifurcation points) the system 
can be in a critical state and can spontaneously transit between two different 
kinds of order (i.e. states). Spontaneous transitions occur only in the critical 
regions Te < T k — A, T k + A > and for L — > oo =>- A — >■ 0, which is usual in 
the critical point. However, it should be noted that in social systems L can 
rarely be treated as infinite and thus A > 0. Within the proposed model, 
'real life' takes place in the critical region. 

Occurrence of spontaneous transitions has been also observed in [14| for 
the Deffuant model with noise, and some analytical results were made, how- 
ever no detailed analysis of this phenomena has been provided. On contrary, 
in this paper an analysis of the distribution of time between transitions (i.e 
state's lifetime) and influence of L p and N on the t(X), have been presented. 
Also influence of the perturbation on ME solutions was investigated. I hope 
it will shed some more light on that BC models analysis. 



5. Appendix 

Derivation of master equation for time evolution of P{0) for model de- 
scribed in this paper is almost the same as in except one small change 
in number of opinions updated in one step. P n {0) is probability density 
function of the opinions at step n and is constructed from the histogram of 
individual opinions O l n . Let's choose two agents i,j to update at step n their 
opinions will be O l n , 0{. Probability that agent i at step n + 1 will adopt 
opinion O is P^ +1 (0) and is given by below formula. 



k + i(o) = f doidoip n (oi)p n (oi)5 to 

+ f dO l n dOiP n {O^P n {Oi)5{0 - 

J\0' n ~0 J n \<T/2 

The independence approximation for the variables O l n , 3 n has been as- 
sumed and it means that P n (O l n , OQ = P n (O l n ) P n (0 3 n ) . To figure out P n+1 (0) 
we have to include interaction between agents with probability ( 1 — p) as well 
as random change of opinion of randomly chosen agent with probability p. 
For Deffuant model equation has form given below. 



O i + O j 



- o l 



(6) 



18 



P n+1 (0) = (1-p) 

+ p 



N - 1 1 
— P.(0) 4- -P.(0) 



(7) 



Now we want to consider this equation for model described in this paper 
for case of L = 1. Because there is only one agent changing it's opinion, the 
part jjPn + i(0) should be neglected and N — 2 should be replaced buy N—l. 



P n+l (0) = (1-p) 
+ P 



N — l 



Pn(0) + -PUl(0) 

N- 1 1 



N 



\ ^(O) + -P a (0) 



(8) 



After replacing P^ +l from Eq. El and simple transformations there goes: 



P n+1 (0) = P n (0) + ^—^[2 [ dO'P n (20 - 0')P n {0') - 

N J\0-0'\<T/2 

P n (0) f dO'P n (0')} + UPJtO) - P n (0% (9) 

J\0-0'\<T 1\ 

For continuum limit P n {0) — > P(0, t) with time t = nSt and 5t = 1/N — > 
as — > oo there is: 



5P(0,t) 
St 



(1 - p) [2 / dO'P(20 - O', t)P(0', t) 

J\0-0'\<T/2 



- P(0,t) dO'P(0',t)]+p[P a (0)-P(0,t)}, (10) 

J\0-0'\<T 

Which is master equation for model described above for L = 1 case, where 
p denotes noise intensity. 
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